Magnetic Excitations in the Spin-Peierls System CuGeOs 



On 
OS 
t— i 

S3 

m 

(N 



i 

CO 

c3 



I 

S3 
O 

o 



> 

so 



O 

Os 
03 



i 

S3 
O 

o 



X 
S3 



G. Bouzerar, A. P. Kampf, and F. Schonfeld 

Institut fur Theoretische Physik, Universitdt zu Koln, 
Ziilpicher Str. 77, Koln 50937, Germany. 

We have calculated the magnetic excitation spectrum in frustrated and dimerized spin 1/2 Heisen- 
berg chains for model parameters which describe the thermodynamics and low frequency spin dy- 
namics in the spin-Peierls system CuGeOz- As a test the chosen model is found to reproduce the 
lowest Raman excitation energy near 30cm -1 in the dimerized phase. We establish the elementary 
triplet and singlet excitation branches below the continuum and estimate the interchain coupling 
effects on the dimerization parameter. 



PACS numbers: 64.70.Kb, 75.10.Jm, 75.50Ee 



Among the current experimental and theoretical ef- 
forts for understanding the magnetic properties of linear 
chain and ladder materials the recent discovery of a spin- 
Peierls (SP) transition at T$p — 1A.3K in the inorganic 
compound CuGeO^ Fj|| has attracted particular atten- 
tion. This interest is partially due to the unique structure 
and available quality of CuGeOj, crystals which allow for 
very detailed studies of the SP phenomenon as well as 
due to the rich phase diagram which evolves in an applied 
magnetic field Q and with doping by substituting Zn for 
Cu H or Si for Ge Q] , respectively. The magnetic prop- 
erties of CuGeO^ arise from Cu 2+ spin 1/2 moments in 
weakly coupled CuOi chains which dimerize below T$p, 
leadingto an alternation of the Cu-Cu distance along the 
chain |g| . This magnetoelastic transition is driven by the 
magnetic energy gain of antifcrromagnetically (AF) cou- 
pled Cu 2+ spins which overcompensates the lattice de- 
formation energy [|| . In the SP ordered phase the copper 
moments form singlet dimers along the chains with an en- 
ergy gap to spin triplet excitations. Experimentally, the 
SP nature of the transition and the spin excitation gap 
have been firmly established by inelastic neutron scatter- 
ing (INS), susceptibility, A-ray, and electron-diffraction 
experiments JtJ. 

The dominant magnetic intrachain coupling between 
Cu 2+ moments in CuGeOz arises from superexchange 
via the bridging oxygens. Although the Cu—O—Cu bond 
angle is near 90° and therefore, by the Goodenough- 
Kanamori rules expected to be ferromagnetic, side group 
effects due to the hybridization of O and Ge orbitals 
lead to an effective AF exchange interaction ||. Fur- 
thermore, Cu — O — O — Cu exchange paths lead to an 
additional sizeable next-nearest neighbor AF exchange 
coupling which frustrates the magnetic interaction be- 
tween the Cu 2+ spins. 

Due to the weak interchain coupling, estimated from 
the dispersion of the magnetic excitations perpendicular 
to the chains to be an order of magnitude smaller than 
the intrachain exchange j^] , a ID approach for modeling 
the magnetic properties of CuGeOs serves as an appro- 
priate starting point. Based on the above arguments we 
use the ID model Hamiltonian 



H=jJ2([l + <5(-l) l ]S 4 • S i+ i + aSi • S 



i+2 



(1) 



where i denotes the sites of a chain with length L and 
Si are 5=1/2 spin operators. J > is the intra- 
chain exchange coupling and a the frustration param- 
eter. The dimerization S accounts for the alternation of 
the Cu — Cu distance along the chain in the dimerized 
phase. For 5 > or for a > a c the singlet groundstate 
of this model is dimerized and a spin gap appears in 
the excitation spectrum. By numerical scaling analysis 
it has been shown that for 5 — the critical frustration 
for spontaneous dimerization is a c as 0.2412 [ fjj5|JIlfl . For 
a < a c and 5 — the system is gapless and renormal- 
izes to the Heisenberg fixed point. For the special case 
2a + 5 = 1 the groundstate is known exactly to be a 
product wavefunction of independent singlet dimers [ fl2"| . 

Considerably different parameter values in Eq. (1) 
have so far been used in previous work on CuGeO^ 
Therefore, as a prerequisite for further CuGeO$ 
specific analysis of the model Hamiltonian, we first fix the 
three parameters by using thermodynamic and INS data. 
With the fixed parameter set we evaluate the Raman in- 
tensity in the dimerized phase and find that the energy 
of the lowest Raman active singlet excitation is smaller 
than twice the triplet excitation gap Aoi . Thus there ex- 
ist clearly at least two elementary excitation branches - 
singlet and triplet, respectively - below the continuum of 
excitations starting at 2 Aoi . Previously it was indicated 
that a singlet excitation might result from a bound state 
of two triplet excitations [p4jjlq| . 

In order to fix the parameters J and a we follow the 
analysis of Biichner et al. |l(| which is based on sus- 
ceptibility (up to lOOOif) Xb(T) data for magnetic fields 
applied along the crystal 6-direction and thermal expan- 
sion data in the uniform (5 = 0) phase; the latter data set 
contains implicit information about the magnetic contri- 
bution to the specific heat. They convincingly conclude 
that clearly a > a c and from the experimental data for 
T > Tsp the exchange couplings are estimated to be 
near J = 160^ and a = 0.35. Indeed, we have veri- 
fied with Lanczos diagonalization techniques applied to 
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FIG. 1. Extrapolated singlet-triplet gap vs. dimerization FIG. 2. Singlet-triplet Aoi and singlet-singlet gap Aoo 

8 for a = 0.35. The dashed line indicates the experimental as a function of 1/L. The symbols are exact Lanczos and the 
results from [HI. dashed lines DMRG results. 



Eq.(l) for 6 = that with these parameters an extraordi- 
narily good fit to Xb(T) in the whole temperature range 
50.fr < T < 1000 AT can be obtained using the g factor 
gt = 2.26 as measured by ESR p7[ . Note that this result 
for a is not in conflict with the work of Castilla et al. |Tl| l 
which favours a < a c . Their arguments were based on 
early neutron scattering data by Harris et al. [Q which 
were recently revised by Martin et al. p9| . 

Given J and a we fix the dimerization parameter 
6 from the requirement that the excitation spectrum 
of Eq.(l) reproduces the low temperature value of the 
singlet-triplet gap Aq^ p w 2.l5meV = 0.156J as mea- 
sured by INS H . We have calculated Aoi on finite chains 
with L < 24 sites as a function of 5 using a Lanczos algo- 
rithm. In order to extrapolate to the infinite chain limit 
L — > oo we use the scaling ansatz 

A 01 (L, S) = A 01 (<5) + - cxp (-J^ffi) (2) 

and an analogous ansatz for the groundstate energy per 
site Eq/L = const. + (B/L 2 ) exp (— £/£). The exponen- 
tial in each case reflects the gapped situation. For the 
groundstate energy the length scale £ is the spin-spin 
correlation length as we have verified by a direct cal- 
culation of the correlation function (S z (r)S z (0)). The 
ansatz Eq.(Q) has been previously motivated and used 
succesfully by Barnes et al. EQ] to study the magnetic 
excitations in gapped spin ladder systems. In Fig.[j] we 
have plotted the extrapolated results of Aoi vs. dimer- 
ization 6 which follow the relation Aoi = c-x + C2<5° with 
the Cross-Fisher exponent a ~ 2/3 as we have verified 
by density matrix renormalization group (DMRG) calcu- 
lations. We observe that the experimental singlet-triplet 
gap value is obtained for 8 w 0.012 for which the ground 
state correlation length is £ ~ 7.3 lattice spacings. 

Having fixed the model parameters we now explore 
their consequence for the Raman intensity. Raman light 



scattering measures the singlet excitations and is thus 
complementary to INS from which the dispersion of the 
triplet excitation has been obtained. Kuroe et al. 
and van Loosdrecht et al. [£lj have observed additional 
peaks in the Raman intensity which appear on cooling be- 
low T$p. The lowest Raman excitation in the dimerized 
phase is observed at 30cto _1 , i.e. slightly below 2 Aoi. 

Within the standard Loudon-Fleury theory of mag- 
netic Raman scattering we write the relevant part of the 
Raman operator for the dimerized phase as 

fffl = A^(S i -S i+1 + 7 S i -S i+2 ) (3) 

i 

for a scattering geometry in which the incoming and scat- 
tered photons are polarized along the chain p2[ ; exper- 
imentally, essentially no scattering is observed in other 
geometries. A is an overall coupling constant while the 
term proportional to 7 arises as a consequence of frus- 
tration [p3| , p2| . The value of 7 depends on microscopic 
details. However, it is expected close to the value of a 
due to their common microscopic origins. With Hr in 
Eq.(||) the scattering intensity is calculated from 

I(uj) ex -I m (0\H R — — 1— - — H R \0) . (4) 
uj + iU + — H + Eq 

We focus our attention on the lowest energy Raman 
excitation in I(u>). In Fig.|| we have plotted Aoi and the 
singlet-singlet gap Aoo as a function of 1/L for our fixed 
parameter set. We present both, exact Lanczos diago- 
nalization results for chains up to 24 sites and DMRG 
data for L < 100. For the DMRG calculations we have 
used the infinite system method with periodic boundary 
conditions keeping 120 states for one target state p4[ ]. 

Fig.^ shows that A o has an unusual nonmonotonic 
scaling behaviour. As a finite size effect for short chains, 
the lowest singlet appears below the triplet excitation; 
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FIG. 3. Scaling of the low excitation energies (in units 
of J) in the k — n momentum subspace. The open (filled) 
symbols are the triplet (singlet) excitations. 



this has also been noted by Riera et al. p5[ . Our DMRG 
and Lanczos data are in excellent agreement for L < 24 
and accurately confirm the ansatz Eq.(|J). Without any 
further scaling fit the DMRG data allow us to read off 
directly the L — + oo value of Aoo ~ 0.232 J. Obviously 
we have Aoo < 2Aoi and Aoo is close to but ~ 13% below 
the experimental result for the lowest Raman excitation 
energy A e Q x p = 30cm" 1 = 0.268J. 

Yet, the strong magnetoelastic coupling as underlined 
by the large magnetostriction observed in CuGeOz |27]] 
implies that magnetic excitation energies will be renor- 
malized due to the coupling to the lattice and a similar 
and presumably larger renormalization is expected from 
interchain coupling. Hence, perfect agreement with the 
experimental data for both the triplet and the singlet ex- 
citation gap can not be expected and we consider the 
above result for Aoo convincing evidence for a consistent 
parameter choice for J, a, and 5. We note that the con- 
trary conclusions of |2q| are not based on a proper finite 
size scaling analysis as in Eq.^|. 

In order to analyze the Lanczos excitation spectrum 
in more detail we consider the strong dimerization case 
6 = 0.2 for a = 0.35 because in this regime the finite 
size effects are very small while we have verified that 
the physics remains essentially unchanged. In fact, the 
ratio R = Aoo/Aoi changes weakly as we increase the 
dimerization parameter S, i.e. R — 1.55 for 5 = 0.2 to 
R = 1.49 for 8 = 0.012. We have plotted in Fig.| the 
scaling of the low energy triplet and singlet excitations 
in the subspace of momentum k = n. We indicate 2Aoi 
in FigJ^ for which we expect to observe the low energy 
edge for the continuum. Surprisingly we observe that also 
the second triplet excitation energy scales below 2Aoi. 

Complementarily we calculate the spectral weight 
of these excitations in the dynamical structure factor 
S(tt,uj) and the Raman intensity I(u>), respectively, as 
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FIG. 4. (a) Scaling of the weight of the 2 lowest Raman 
singlet excitations calculated with 7 = 0.3 in Eq.^. (b) Spec- 
tral weight of the low energy triplet excitations in S(tv,u) vs. 
1/L. For the lowest excitation the weight is multiplied by 0.2. 
In (a) and (b) the same symbols are used as in Fig|| 

given by Wi — |(i|O|0)| 2 where |0) is the groundstate 
and \i) is an excited state with momentum k = n. 
The operator O is Hr/Vl for the singlet and S z (k) = 
{\/\/T,)^2 l exp{ikl)Sf for the triplet excitations. In 
Fig.||b we show results for Wi of the 3 lowest triplet ex- 
citations as a function of 1/L. It appears clearly that Wi 
for the two lowest triplets scales to a finite value, identify- 
ing them as elementary excitations below the continuum. 
Contrary, Wi for the third excitation in the triplet sub- 
space scales to while its energy extrapolates to 2Aoi, 
i.e. the onset of the continuum. This is also reflected 
by the kink in the scaling data for L = 16. In Fig^a we 
perform the same scaling of Wi for the singlet excitations 
in the Raman intensity evaluated for 7 = 0.3 in Eq. ([?]). 
The conclusions remain unchanged if 7 is varied. Wi of 
the lowest singlet excitation obviously scales to a finite 
value while the weight of the next higher energy singlet 
excitation scales to 0. 

From this finite size scaling analysis we conclude that 
there are 3 elementary excitations, 2 triplet and 1 sin- 
glet, for S = 0.2; the two lower excitations clearly persist 
when 5 is decreased to 0.012 - the model parameter value 
relevant for CuGeO-s . The singlet excitation at k = (or 
k — ir) is identified with the 30cm -1 peak in the Raman 
intensity while the lower triplet corresponds to the spin 
excitations measured in INS. For the existence of the sec- 
ond triplet near the edge of the continuum, however, a 
definite conclusion could not be reached for <5 = 0.012. 
The results for 5 = 0.2 are summarized in Fig.|| which 
shows the dispersions for the 3 elementary excitations. 
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FIG. 5. Dispersion ui(k) (in units of J) of the elementary 
triplet and singlet excitations. The continuum from all spin 
sectors above 2 Aoi (fc) is indicated by the shaded area. 

So far we have used a strictly ID approach to describe 
the magnetic excitations in the SP phase of CuGeO^. 
Yet, for a more realistic model the interchain coupling via 
the oxygen atoms that two adjacent CuO^ chains have in 
common cannot be ignored. Indeed, the 6-direction dis- 
persion of the triplet excitation measured in INS is not 
small and its energy varies from its maximum value of 
5.8meV" to 2.1meV at the Brillouin zone (BZ) center 
In a first attempt to include interchain coupling effects 
into an effective ID spin Hamiltonian we may determine 
the dimerization parameter from the BZ averaged triplet 
excitation gap w 3.95meV rather than from its BZ cen- 
ter value used above. Keeping J and a fixed as before, 
because they are determined from data in the uniform 
phase, then leads to an estimate S ~ 0.036, i.e. an in- 
crease by a factor of 3. Assuming a similar 6-direction 
dispersion also for the singlet excitation branch we find 
that the BZ center value of the singlet gap remains al- 
most unchanged and thus our conclusions are still valid. 

In conclusion, we have fixed the parameters of a frus- 
trated and dimerized Heisenberg chain model to describe 
the thermodynamics and the low energy magnetic exci- 
tations at low temperatures in CuGeOj,. The internal 
consistency of the parameter choice is verified by almost 
quantitatively reproducing the energy of the lowest sin- 
glet excitation in the dimerized phase as measured by 
Raman light scattering. Furthermore, we have shown 
the existence of at least 1 elementary triplet and 1 single 
excitation branch below the continuum. While the effects 
of weak interchain coupling on the spin dynamics still re- 
main to be explored the J—a — S model Hamiltonian with 
fixed parameters will serve as a basis for studying doping 
and magnetic field dependent phenomena in CuGeO^. 
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